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Abstract 

The critical sector of strong interactions at high temperatures is explored in the 
frame of two complementary paradigms: Statistical Bootstrap for the hadronic phase 
and Lattice QCD for the Quark-Gluon partition function. A region of thermodynamic 
instability of hadronic matter was found, as a direct prediction of Statistical Bootstrap. 
As a result, critical end point solutions for nonzero chemical potential were traced in 
the phase diagram of strongly interacting matter, consistently with chiral QCD and 
recent lattice calculations. The relevance of these solutions for current and future 
experiments with nuclei at high energies is also discussed. 
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1. Introduction 

Quantum Chromodynamics is unquestionably the fundamental theory of strong inter- 
action. However, the nonperturbative aspects of QCD belong still to a field of intense 
investigation. In fact, most of the properties of the hadronic world cannot be extracted 
yet from first QCD principles and the partition function of interacting hadrons produced 
and thermalized in high-energy collisions can only be determined within specific models of 
varying degrees of limitation. 

In this paper we employ Statistical Bootstrap, including volume corrections for the finite 
size of hadrons, in order to describe the hadronic phase. The virtue of this description 
is associated with the fact that, among Statistical models of hadrons, it is the only one 
which predicts the onset of a phase transition in strongly interacting matter and therefore 
it is compatible with the properties of the QCD vacuum at high temperatures [1]. The 
aim of this work is to pursue the search for this compatibility with QCD even further, 
asking whether Statistical Bootstrap of hadrons is compatible with the existence of a critical 
endpoint in strongly interacting matter, at high temperature and finite (nonzero) baryonic 
chemical potential. The existence of such a critical point in the phase diagram is required 
by QCD in extreme conditions [2] and it is the remnant of chiral QCD phase transition [2]. 

The basic ingredients in our approach are (a) the hadronic partition function extracted 
from the equations of Statistical Bootstrap and (b) the equation of state of the quark-gluon 
phase given by recent QCD studies on the lattice [3-6] . Our principal aim from the matching 
of these two descriptions is to trace the formation of a critical point in the quark-hadron 
phase transition with a mechanism compatible both with Statistical Bootstrap and Lattice 
QCD. 

In Section 2 the principles and basic equations of Statistical Bootstrap, leading to the 
partition function of hadronic phase, are briefly summarized and the appearance, within 
this framework, of a thermodynamic instability is discussed in detail. This instability is the 
origin of the formation of a critical point in the bootstrap matter. In Section 3 the partition 
function of the quark-gluon phase is extracted from recent lattice calculations of the pressure 
of QCD matter [3]. In Section 4 the above two descriptions of strongly interacting matter 
are exploited in a search for a critical point in the phase diagram, the location of which 
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is fixed by a set of equations (21-25). Finally, in Section 5 our conclusions including the 
limitations of our approach are presented whereas in the Appendix certain technical points 
are presented concerning the evaluation of the quark-gluon partition function (part A) and 
the full form that acquires the equation of maximum hadronic pressure (part B). 

2. Thermodynamic Instability in Statistical Bootstrap 

The main attribute of the Statistical Bootstrap (SB) [7-10] is that it describes thermo- 
dynamically hadronic systems which are interacting. Generally as far as its thermodynamic 
behaviour a system of strongly interacting entities may be considered as a collection of par- 
ticles with the explicit form of the complex interaction among themselves. Another way 
to deal with the problem is to assume that the strongly interacting entities form a greater 
entity with certain mass. Then, if this mass is known, the greater entities may be treated 
as non-interacting and the simple thermodynamic description of an ideal system of parti- 
cles may be applied. The SB adopts the second approach and assumes that the strongly 
interacting hadrons form greater clusters, called "fireballs" . The problem then is moved to 
determine the mass of these clusters or to evaluate their mass spectrum f(m 2 )dm 2 , which 
is equal to the number of discrete hadronic states in the mass interval {m,m + dm}. The 
solution is based on the assumption that the fireballs may consist of smaller fireballs with 
the same mass spectrum or "input" particles which may not be divided further. The integral 
bootstrap equation then reads [11-14] 

00 1 r ( n \ 
B(p 2 )f(p 2 , b,s,...)= g b , s ,.B{p 2 )W ~ ml,...) + E - M> " EH • 

> . ' n=2 Hl J \ i=l ) 

input term 

(n \ / n \ n 

*-EMEM a -5>---n b(pMpI h, 8i , . . .)d% . (i) 
1=1/ { Sl } \ i=l ) i=l 

Equation (1) also imposes conservation of four- momentum p and any kind of existing 
additive quantum number (baryon number b, strangeness s, etc.) between a fireball and its 
constituents fireballs. The masses m^... are the "input" particle masses which constitute 
the smaller fireballs that can be formed. 

The appropriate volume of a fireball is considered to be carried with it, = The 
imposition of the requirement that the sum of volumes of the constituent fireballs has to be 
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equal to the volume of the large fireball, as well as, momentum conservation lead to the fact 
that all fireballs posses the same volume to mass ratio. This ratio can be connected to the 
MIT bag constant, ^ = j^. The correct counting of states of a fireball involves, apart from 
its mass spectrum which is of dynamical origin, a kinematical term B 

R(J\ - 2WP » - Wm (2) 
B(P } " W " (2^ ■ (2) 

It can be noted in (1) that the mass spectrum is accompanied by the term B. Thus the 
bootstrap equation remains unchanged if f and B are redefined in such way so 

~ A — ry A — /y 
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where a is a free real number and m a free parameter with the same dimensions as the mass. 
The effect of such transformations will become apparent when we turn to the thermodynamic 
description of the system of fireballs. 

The bootstrap equation can be simplified after performing a series of Laplace transfor- 
mations to acquire the form 

¥>(/?, {A}) = 2G((3, {A}) - exp[C(/3, {A})] + 1 . (4) 

In the last equation G(f3, {A}) is the Laplace transform of the mass spectrum with the 
accompanying kinematical term 
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and ip(f3, {A}) the Laplace transform of the input term 
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The masses mb qs \ s \ correspond to the masses of the input particles, which in this paper will 
be all the known hadrons with masses up to 2400 MeVs, the gb qs \ s \ are degeneracy factors due 
to spin and the A's are the fugacities of the input particles. Here we have used the extended 
version of SB [14], where the states are characterised by the set of fugacities relevant to 



baryon number, A&, electric charge, X q , strangeness, A s and partial strangeness equilibrium, 
7 S (in the following we shall refer to this set of fugacities as {A}, for short). 

The bootstrap equation defines the boundaries of the hadronic phase since it exhibits a 
singularity at the point 

^(/3,{A}) = ln4-l. (7) 

From the physical point of view the singularity is connected with the behaviour of the 
mass spectrum as the mass tends to infinity 

f (m\ {A }) 2C({A})m-«- 1 exp[m/T ({A})] , (8) 

where (3 = T~ l and [3* corresponds to the maximum inverse temperature. After a certain 
point, as temperature rises, it is more preferable for the system to use the given energy in 
producing more hadronic states (since their number rises exponentially) than in increasing 
the kinetic energy of the already existing states. 

In order to turn to the thermodynamics it is necessary to consider the fireball states in an 
external volume V ext . This volume must be distinguished from the physical volume which is 
carried by each fireball. The partition function of the pointlike interacting hadrons is then 

\nZ(V ext ,f3,{\}) = 

E *r £ X t E As E 71 s ' 1 / -^^f\p 2 ,b>, q >,s>,\s>\)e-^dp*. (9) 

b'=-oo q'=-oo s'=-oo \s'\=0 \ Zn ) 

Every choice for the mass spectrum (or for the exponent a) in (3) leads to different partition 
function and so to different physical behaviour of the system. The usual SB choice was 
a — 2, but more advantageous is the choice a — 4. With this choice a better physical 
behaviour is achieved as the system approaches the hadronic boundaries. Quantities like 
pressure, baryon density and energy density, even for point-like particles, no longer tend to 
infinity, as the system tends to the bootstrap singularity. It also allows for the bootstrap 
singularity to be reached in the thermodynamic limit [15], a necessity imposed by the Lee- 
Yang theory. Another point in favour of the choice a = 4 comes from the extension of SB to 
include strangeness [11,12]. The strange chemical potential equals zero in the quark-gluon 
phase. With this particular choice of a, /i s acquires smaller positive values as the hadronic 



4 



boundaries are approached. With the choice a = 4 the partition function can be written 
down and for point-like particles it assumes the form 

ABV 



A Hi/ poo 

In Z(V ext , (3, {A}) = — J x'G(x, {\})dx = Vf(P, {A}) 
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For a = 4 the input term acquires the form 
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where "a" represents a particular hadronic family, u i" the hadrons in the family with different 
masses and B is the energy density of the vacuum (MIT bag constant). 

By including corrections due to the finite size of hadrons (Van der Waals volume correc- 
tions) the repulsive part of the interaction is taken into account. The negative contributions 
to the volume can be avoided if the following grand canonical pressure partition function is 
used 



tt(£,/3,{A}) 



(12) 



£ - f(P + {A}) ' 

which is the Laplace transformed partition function with respect to volume [16] (£ is the 
Laplace conjugate variable of the volume). All values of £ are allowed if Gaussian regular- 
ization is performed [16]. In the following we shall consider systems without any external 
forces acting on them and so we shall fix £ = [15,16]. The density and the pressure P of 
the thermodynamic system can be obtained through the pressure grand canonical partition 
function (12) for £ = 



^ G (£ = 0,/3,{A}) = A 
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AB 



dp 
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where A is the fugacity corresponding to the particular density, and 
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Though volume is no longer an active variable of the system it can be calculated for 
given baryon density and v B (evaluated through (13)) and baryon number < B > which is 
a conserved quantity. The volume would be retrieved through the relation 

< B > 



< V >-- 



(15) 



With the use of SB in order to describe interacting hadronic systems we can trace the 
possibility of a phase transition. The study of the pressure-volume isotherm curve is then 
necessary. When this curve is calculated a region of instability is revealed. In fact, this 
curve has a part (near the boundaries of the hadronic domain) where pressure decreases 
while volume decreases also (see Fig. 1). Such a behaviour is a signal of a first order phase 
transition which in turn can be mended with the use of a Maxwell construction. 

This behaviour is due to the formation of bigger and bigger clusters as the system tends 
to its boundaries in the phase diagram. In that way the effective number of particles is 
reduced, resulting, thus, to a decrease of pressure. To show that this instability in the 
P — V curve is the result of the attractive part of the interaction included in the SB we 
shall calculate a similar curve using the Ideal Hadron Gas (IHG) model with Van der Waals 
volume corrections (repulsive part of interaction). The logarithm of the partition function 
of IHG (corresponding to (10)) is 

In Z p IHG (V, (3, {A}) = Vf p IHG (P, {A}) = 

E[ A a({A}) + A a ({A})^]5: fe m ai K 2 (/3m a ,) , (16) 

where are degeneracy factors due to spin and isospin and the index a runs to all hadronic 
families. This function can be used in eq. (12) to calculate the Ideal Hadron Gas (IHG) 
pressure partition function in order to include Van der Waals volume corrections. The result 
is that the pressure is always found to increase as volume decreases, for constant temperature, 
exhibiting no region of instability and so no possibility of a phase transition. 

The comparison of SB with the IHG (with volume corrections) is displayed in Fig. 1, 
where is is the normal nuclear density is = 0.14 fm~ 3 . In both cases (SB or IHG) the 
constraints < S >= (zero strangeness) and < B >= 2 < Q > (isospin symmetric system, 
i.e. the net number of u and d quarks are equal) have been imposed. Also strangeness is 
fully equilibrated which accounts to setting j s — 1. 

3. The partition function of Quark Matter 

Having a description of the hadronic phase at hand, it is necessary to proceed with the 
thermodynamical behaviour of the quark-gluon phase. Lattice calculations of the pressure 
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of the quark-gluon state have been performed at finite chemical potential in [3,5]. These 
publications include calculations for rather heavy u, d quark masses. The mass of the u, d 
quarks is about 65 MeV and the strange s quark 135 MeV [3,5]. The calculated pressure of 
the quark-gluon phase (P/T 4 ) at \i# = is plotted against the ratio of temperature to the 
transition temperature of quark matter at zero baryon chemical potential T/T c in Fig. 2 of 
[3]. The temperature T c will be denoted as T qgp in the following. The results of this graph 
are extrapolated to the continuum limit by multiplying the raw lattice results with a factor 
c p = 0.518 [3]. 

The lattice calculations for finite chemical potential can be summarised in Fig. 3 of 
[3], where the difference of pressure at non-zero chemical potential and the pressure at zero 
chemical potential {Ap/T 4 = [P{fi ^ 0,T) - P(/i = 0,T)]/T 4 ) is plotted against T/T c . 
Again the results of this graph are extrapolated to the continuum limit by multiplying the 
raw lattice results with a factor = 0.446 [3]. In this graph there are plotted five curves 
which correspond to baryon chemical potential of 100, 210, 330, 410, 530 MeV. 

With the use of Figs. 2, 3 in [3], it is possible to calculate in principle the pressure of 
the quark-gluon phase at any temperature and baryon chemical potential. The pressure is 
important, because knowledge of the pressure is equivalent to the knowledge of the partition 
function of the system in the grand canonical ensemble 

\nZ QGP (V,T,^ B ) = ^P(T^ B ) (17) 

In order to have a complete description of the dependence of the pressure on the temper- 
ature and the chemical potential we use two sets of fitting functions. For constant chemical 
potential the pressure as a function of T/T c is fitted through 

/<*> = ; ; ,f, + ; * ^ («0 



X L 
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where a,i,bi,Ci,di, fi (i = 1, 2) depend on while for constant temperature the correspond- 
ing fit of the pressure as a function of /ib is given by 

g{x) = a + bexp(cx d ), (19) 

where a, b, c, d depend on the temperature ratio T/T c . The fitting procedure has to be 
performed in a self-consistent way and subsequently it is straightforward to evaluate the 



partition function as well as its derivatives with respect to /j,b and T at any given point 
(Ti, /ib i)- In particular to evaluate physical observables connected with the partition func- 
tion and drive numerical routines the partial derivatives of the pressure up to second order 
with respect to temperature and fugacity have to be evaluated. These derivatives are then 
given in part A of the Appendix. 

In Fig. 2 we have reproduced the quark-gluon pressure as a function of the temperature 
for constant baryon chemical potential. The squares are points directly measured from the 
graphs of Fodor et. al and the lines represent the calculation with the fits which has been 
performed on these points, via eq. (18). Fig. 3 is a graph similar with Fig. 2, but we have 
focused on the area which is useful for our calculations, the area where the matching with 
the hadronic phase will be performed. Fig. 4 is a reproduction of the Fodor et. al quark- 
gluon pressure as a function of the baryon chemical potential for constant temperature. The 
necessary fits have been performed with the use of eq. (19). 

4. The critical point in the phase diagram 

After developing the necessary tools to handle the thermodynamic description for the 
quark-gluon phase as it is produced from the lattice, we can search for the possibility of a 
quark-hadron phase transition. We shall choose £ = [15] and so the only parameter left, 
in relevance with the hadronic side, would be the maximum temperature at zero baryon 
chemical potential T hg- We can, also, allow the critical temperature of quark-gluon state 
at zero baryon chemical potential Tqqgp as free parameter. The fact that To hg and Tqqgp 
are treated as two separate parameters does not include a contradiction, since T hg is the 
maximum temperature that the hadronic phase can exist, in the bootstrap formalism, at 
fiB = 0. This does not exclude the fact that the phase transition to the quark-gluon phase 
can take place at a smaller temperature at zero chemical potential. So it can happen that in 
general Tqqgp < Tq hg and this is the case as it will become evident in the following. Also, 
the region of small baryon chemical potentials belongs, according to the predictions of QCD, 
to the cross over region, where no distinct separation between the hadronic and the quark 
state exists. 

Then, if values for these parameters are chosen, it is possible to calculate for a specific 
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temperature the pressure isotherms of Hadron Gas and QGP. Assuming that the baryon 
number is a conserved quantity to both phases, the equality of volumes is equivalent to the 
equality of baryon densities, and so the connection of the isotherms of the two phases is 
possible through the relation 

< V HG >=< V QGP >^ ^— = ^— ^u BHG = v B QGP (20) 

Vb HG v B QGP 

The graph of the pressure-volume isotherm can be drawn if the pressure is plotted against 
the inverse baryon density. Then at the point where the isotherms of the two phases meet 
we have equal volumes for equal pressures. 

Tracing the point where the isotherms of two phases meet, what is found is that at a 
low temperature the QGP and SB Pressure- Volume isotherms intersect at a point where 
the Hadron Gas pressure is decreasing while volume decreases. The resulting pressure- 
volume curve includes an unstable part which has to be repaired through a suitable Maxwell 
construction. This curve includes a region where a first-order transition takes place. As 
temperature increases, a certain temperature is found for which the QGP and SB isotherms 
meet at a point which is, also, the maximum Hadron Gas pressure for this temperature. In 
that case no Maxwell construction is needed and since this point is located at finite volume 
or not zero baryon density (equivalently not zero chemical potential) it can be associated 
with the QCD critical point. As temperature rises more, the QGP and SB isotherms meet at 
a point which corresponds to even greater volume. Again no Maxwell construction is needed 
and since the resulting pressure-volume isotherm always increases while volume decreases 
the meeting point of the two phases belongs to the crossover region. 

A graph that summarises the situations met in the pressure volume isotherms of hadronic 
and quark systems is Fig. 5. In this figure the hadronic isotherms have been calculated for 
T jjg — 177 MeV (MIT bag constant 5 1/4 ~ 222 MeV) while the quark-gluon isotherms for 
T qgp ~ 152 MeV. For T = 161.3 MeV a Maxwell construction is needed to remove the 
instability from the resulting curve. The horizontal line defines a partition with two equal 
surfaces (shaded) and represents the final pressure- volume curve after the completion of the 
Maxwell construction. At the temperature T = 162.1 MeV the two curves meet at the point 
of maximum hadronic pressure and so a critical point is formed at finite volume. At a greater 
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temperature (T ~ 162.3 MeV) the quark pressure-volume curve meets the hadronic one at 
a point away from the area of instability (crossover area). 

To locate the critical point numerically with the use of the lattice partition function, for 
given parameters T HG and T q GP , the conditions have to be determined for which the SB 
pressure is equal to the QGP pressure at the same volume, corresponding to the maximum 
SB pressure. Setting the factor of partial strangeness equilibrium 7 S = 1 a hadronic state 
is characterised by the set of thermodynamic variables (T, A M , A^, A s ), while a quark-gluon 
state evaluated by the lattice [3] is characterised by the two variables (T, A'). The u and 
d quarks are characterised by the same fugacity A^ = A^ = X' q = X' B 3 . To evaluate the 
unknown variables we have to solve the following set of non-linear equations. 

v B sb(T, X u , X d , A s ) = v B QG p(T, X' q ) (21) 

P SB (T, X u , X d , X s ) = Pqgp(T, X' q ) (22) 

Psb(T, X u , Xd, A s ) = P S b(T, X u , Xd, X s ) max (23) 

(S) SB (T,X u ,X d ,X s )=0 (24) 

( B ) SB (T, X u , X d , Xs) - 2 ( Q ) SB (T, X u , X d , X s ) = (25) 

Eq. (21) accounts for the equality of the baryon densities of the two phases which is equivalent 
to the equality of volumes, since the baryon number is a conserved quantity. Eq. (22) is the 
equality of pressures of Hadron Gas and QGP. Eqs. (21), (22) determine the point where 
the two pressure curves meet. Eq. (23) requires that the meeting point of the two phases for 
a certain temperature is equal to the point which maximizes the Hadron Gas pressure and 
so this meeting point is the critical point. The form of this equation is discussed in detail 
in the Appendix. Eq. (24) imposes strangeness neutrality in the hadronic phase. Eq. (25) 
imposes isospin symmetry to the hadronic system. The solution for systems without such 
an attribute is available for the hadronic side but it is of no use in the current study since 
the lattice results with only one fugacity for the quark system have included the assumption 
of isospin symmetry. 

Before presenting results on the existence and location of the critical point, certain limi- 
tations on the available parameters T hg and Tqqgp are in order. To have a critical point 
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which could be associated with QCD it is necessary for this point to be located at non zero 
chemical potential. So to impose limits on the values of these parameters we have to eval- 
uate for which choices the critical point corresponds to zero chemical potential. At hb — 
the constraints (S) SB = and (B) SB — 2 (Q) SB = may be solved analytically to give 
\ u = Xd = A s = 1. Then, v B SB = and consequently vbqgp has to be set to zero, leading 
to X' q — 1. Only, the temperature T, is left undetermined from the set of thermodynamic 
variables of the system. For given parameter T hg one can, then, solve the set of eqs. (22) 
and (23) to determine the value of the set (T;T q GP ). The results of this calculation are 
depicted on Fig. 6, where the solid line represents the connection of the two parameters 
for a critical point at /ib — 0. It is easy then to check that the region on the upper left 
of this curve represents choices of the parameters that drive the critical point on positive 
baryon chemical potential, whereas the region on the bottom and right represents choices 
of the parameters for which the set of eqs. (21)-(25) has no solution. It is evident from 
the slashed line, that it is not possible to have solution for the critical point for the choice 
To hg — Tqqgp- The difference, though, DT = T HG — T QGP between the two parameters 
has a minimum, which is about 5.1 MeV and it occurs for the choice T 0HG = 140.6 MeV 
(B 1 / 4 = 136 MeV) and T 0QGP = 135.5 MeV. 

The solutions for the position of the critical point in the (T, /j, B ) plane for specific choice of 
the parameters are presented in Fig. 7. The solid lines represent solutions with fixed T q GP 
and varying T 0HG , whereas the slashed lines represent the opposite situation. It is evident 
that an increase of T HG at constant T q GP leads the critical point to higher temperature 
and higher baryon chemical potential. An increase of T q GP at constant T HG drives the 
critical point at higher temperatures but at smaller baryon chemical potential. 

Recent lattice QCD studies offer, apart from the quark-gluon pressure which has been 
a basic ingredient in our approach, important results on the existence and location of the 
critical point itself. In [17] the critical point is found to reside at T cr . p . = 160 ± 3.5 MeV and 
Hb = 725 ± 35 MeV, with T q GP = 172 ± 3 MeV. These calculations have the drawback 
that have been performed with rather unphysical mass for u and d quarks which has a value 
about four times the physical value. Improved calculations have been performed in [18], 
where the light quark masses have decreased by a factor of 3 down their physical values. 
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The critical point is found now (with T QGP = 164 ± 3 MeV) to be at T cr . p . = 162 ± 2 MeV 
and fi B = 360 ±40 MeV, a value which is considerably reduced with respect to the previous 
one. This point is depicted on Fig. 7 with the full circle. In our study it can be approached 
with the choice of the parameters T HG = 152.5 MeV and T q GP = 177.5 MeV. 

In Fig. 8 we present our solution for the QGP-hadron transition line in the (T, fx B ) plane. 
The parameters are chosen so as to make the critical point to be in the same location as 
the lattice solution. The circle represents our solution for the critical point and the star the 
lattice result [18]. The solid thick line represents the boundaries of the hadronic phase (as set 
from the bootstrap singularity) which is close to the first order critical line. The solid thin 
line is the first order transition line evaluated in [17] with large u,d quark masses. The dotted 
line represents the crossover in this calculation. The corresponding solution for the critical 
point is also depicted, at a large value of baryon chemical potential. For completeness, 
early solutions for the critical point [2] based on different approximate theories (Nambu- 
Jona-Lasinio (NJL) model or a random matrix (RM) approach) are also shown in the same 
figure (NJL-RM). We have, also depicted on the same graph the freeze-out points of different 
experiments. It is evident that the recent lattice calculations [18] set the critical point to a 
location easily accessible by experiments. 

5. Conclusions 

Statistical Bootstrap presents a more accurate description of the hadronic phase than the 
ideal Hadron Gas, since it includes in a self consistent way the interaction among hadrons. 
This interaction is crucial to investigate critical phenomena in connection with the state 
of quark-gluon plasma. Among the predictions of the bootstrap model is the limitation of 
the hadronic phase and the forming of an instability in the pressure-volume isotherm near 
the hadronic boundaries. This instability can be connected with a first order quark-hadron 
phase transition and the existence of a critical end point in the strongly interacting matter. 

The accurate partition function of the quark-gluon phase is available from lattice calcu- 
lations, which include, though, large values of the light quark masses. From these results 
the lattice partition function of the quark state, as well as, all the necessary derivatives can 
be calculated, allowing the evaluation of any physical observable. 
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The joining of the SB and the lattice partition function for the haronic and the quark 
state respectively, allows for the determination of a critical point at finite baryon chemical 
potential which can be related to the critical point of QCD. 

More recent lattice calculations [18] drive the position of the critical point to smaller 
values of baryon chemical potential as the values of u — d quark masses approach their 
physical values. It is interesting that the current location is situated in the (T, hb) plane in 
a region easily accessible by the freeze-out conditions of experiments at the CERN/SPS. 

With suitable choice of the temperature at zero density (T ) for both the hadronic and 
the quark-gluon system the position of the critical point, determined in our study through 
the statistical bootstrap and the lattice pressure, can coincide with the one directly found 
from lattice calculations. 

In a previous work [19,20] a similar solution was found with the use of a simplified 
partition function for the quark-gluon system, based on the MIT bag model. Therefore, 
the basic mechanism in our approach for the formation of a critical point in the strongly 
interacting matter is not associated with the details of the quark-gluon partition function 
but mainly with the instability of hadronic matter revealed by Statistical Bootstrap. In 
particular the local maximum of pressure in the P — V diagram (fig. 5) of hot hadronic 
matter lies in the origin of the formation mechanism of the critical end point. 

In the absence of a final theory of strongly interacting matter, providing a unified par- 
tition function for both phases, our approach is necessarily an approximate treatment with 
drawbacks and limitations. In fact the continuity of the derivatives of P — V isotherms in the 
crossover area (T > T c ) is not guaranteed (fig. 5) and the power law at the critical tempera- 
ture (T = T c ), dictated by the isotherm critical exponent S, is not revealed in this treatment. 
Despite these limitations, the overall consistency of the critical sector of strong interactions 
at high temperatures, emphasized in this approach, provides us with extra confidence for 
the existence and the location of the QCD critical end point and makes its experimental 
discovery even more challenging [21]. 
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Appendix 

A. The fitting procedures on curves of constant temperature and chemical potential 
allows us to evaluate the derivatives with respect to T for constant hb and with respect to 
Hb for constant T. Physical observables, however, are given as derivatives with respect to 
temperature for constant fugacity or with respect to \ B for constant T. The evalution of 
the latter (for the pressure) is given by 
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OX 
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dp 



dP 
df 



As 
d 2 P 



dP 
df 



T dfiB 
dP 
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dX f 
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Xb^b 

J^2 



(Al) 



dP 
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dP 
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d 2 P 2/i B d 2 P 



d\i B &T T 
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d\ B dT dTdfiB X B dfi 
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dfi% 
dP 



T 



X b dfi B 



/4 

J^2 ' 



Xb 



(A4) 



(A5) 



In eqs. (A. 4), (A. 5) where the 2nd partial derivative with respect to two different variables 
of P appears, the pressure is considered as a function of these two variables. 

B. The requirement of eq. (23) has to be fulfilled for a certain temperature T = T\ and in 
the presence of two constraints g\ = ( S ) SB (7\, A„, X d , A s ) = and g 2 = ( B ) SB (T±, X u , Xd, X s ) — 
2 ( Q ) SB (Ti, X u , Xd, A s ) = 0. The temperature 7\, since the maximum of pressure is found 
for a certain isotherm, may not be considered in the following as an active variable. 



dPsMTu X u , Xd, A s ) 



dP SB = 



dPsB,, . dP S B, x . dP S B, x 
dX u + — — dX d + ^rr — dX s 



d\. 



dPsB OPsBdXd dP SB dX, 
dX u dXd dX u dX s rfA, 



dX d 
■ = 



<9A, 



[B.1) 



As far the constraint g\ is concerned, we have 



gi{T u X u ,X d ,X s ) = d 9l = 

dgi dX d 



dg 



dX u 
dgi dX s 



-dX,, 



T7^dX d + 77T-«A S = 

dXd oX s 



dX d dX u dX s dX u 



dgi 
dX u 



(B.2) 
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Similarly for the constraint g 2 we have 

dg 2 d\ d dg 2 d\ s dg 2 



(B.3) 



dX d d\ u dX s d\ u dX u 

Eqs. (B.2) and (B.3) may considered as a system of two equations which can be solved to 
determine d\d/d\ u and d\ s /d\ u 

d\ d 1 / dgi dg 2 dg 2 dg 1 \ 



d\ u D \ dX s dX, dX s <9A„ 



(BA) 
(B.5) 



d\ a _ 1 / dg x dg 2 dg 2 dg x \ 
d\ u D \dX u d\d d\ u d\ d j' 

with 

D = dg^dg^ _ dg^dg^ 

dx d dx s dx d dx s - 1 • ' 

Eqs. (B.4) and (B.5) may now be inserted to eq. (B.l) to give 

dP S B dP S B ( dgi_dg2_ _ dg 2 dg ± \ j_ dP SB ( dg^dg-^ _ dg^dgA j_ = 
dX u + dX d \dX s dX u dX s dX u ) D + dX s \dX u dX d dX u dX d ) D 1 ' ' 

Eq. (B.7) is the form of eq. (23) for the maximum hadron pressure in a certain isotherm. 
The main contribution comes from the term dPsB/9X u and the result for the maximum 
pressure is near to the maximum evaluated with the full equation (B.7), which explains why 
only the first term was used in [19,20]. 
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Figure Captions 

Fig. 1 Isotherm pressure-volume curve for SB and IHG (both with Van der Waals volume 
corrections using the pressure ensemble). The SB curve is drawn for T HG = 173 MeV 
(5 1 / 4 = 210 MeV). 

Fig. 2 The pressure of the quark-gluon state divided by T 4 versus the ratio T/T 0G . GP , for 
constant baryon chemical potential. The lines from bottom to top correspond to gradually 
increasing values of The squares represent direct measurement from the figs 2, 3 in [3] 
which depict the lattice calculation and the lines our fits on these points. 

Fig. 3 Similar graph with Fig. 2. The pressure of the quark-gluon state is divided by Tqq GP 
and the graph is focusing on the region where the matching with the hadronic state will take 
place. 

Fig. 4 The pressure of the quark-gluon state divided by Tqq GP versus the baryon chemical 
potential, for constant values of the ratio T/T q GP . The squares represent direct measure- 
ment from the figs 2, 3 in [3] which depict the lattice calculation and the lines our fits on 
these points. 

Fig. 5 Three isotherm pressure-volume curves for Hadron Gas (using SB) and QGP phase 
(using the lattice pressure of [3]). The low temperature isotherm needs Maxwell construction, 
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the middle temperature isotherm develops a critical point and the high temperature isotherm 
corresponds to crossover. 

Fig. 6 Connection between the parameters T hg and T q GP (solid line) in order for the 
critical point to be developed at zero baryon chemical potential. The region which leads 
to critical point at positive chemical potential (upper left) and the region with no solution 
(bottom right) are depicted. 

Fig. 7 Position of the critical point at the (T,/ib) plane for constant values of T hg and 
varying T q GP (solid lines) and constant values of T q GP and varying T 0HG (slashed lines). 
The lattice calculated critical point in [18] is displayed by the full circle and the cross. 

Fig. 8 Position of the critical point (open circle) (at the (T, /jl b ) plane for the specific choice 
of parameters T hg = 177.5 MeV and T q GP = 152.5 MeV as well as the first order part 
quark-hadron transition (bootstrap singularity-thick line). The stars represent the lattice 
calculated critical points in [17] and [18]. The thin line is reproduced from [17] and represents 
lattice calculation for the transition line. There are, also, depicted freeze-out points from 
different experiments. 
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Fig. 1 
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